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^ Abstract 

Q The two-body break up of the deuteron is studied at high Q kinematics, with main motivation 

to probe the deuteron at small intcrnuclcon distances. Such studies are associated with the probing 
^ of high momentum component of the deuteron wave function. For this, two main theoretical issues 

Q have been addressed such as electromagnetic interaction of the virtual photon with the bound 

nucleon and the strong interaction of produced baryons in the final state of the break-up reaction. 
Within virtual nucleon approximation we developed a new prescription to account for the bound 
' (— ( ' nucleon effects in electromagnetic interaction. The final state interaction at high kinematics is 

calculated within generalized eikonal approximation (GEA). We studied the uncertainties involved 
Q in the calculation and performed comparisons with the first experimental data on deuteron elec- 

^ trodisintegration at large Q^. We demonstrate that the experimental data confirm GEA's early 

I— —I prediction that the rescattering is maximal at ~ 70*^ of recoil nucleon production relative to the 

^ momentum of the virtual photon. Comparisons also show that the forward recoil nucleon angles 

^ are best suited for studies of the electromagnetic interaction of bound nucleons and the high mo- 

mentum structure of the deuteron. Backward recoil angle kinematics show sizable effects due to 
the A-isobar contribution. The latter indicates the importance of further development of GEA to 
^ account for the inelastic transitions in the intermediate state of the electrodisintegration reactions. 
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I. INTRODUCTION 



Two-body electro-disintegration of the deuteron at high represents a powerful tool 
for studying one of the most fundamental issues of nuclear physics such as nuclear forces at 
intermediate to short distances. Despite all the successes in constructing interaction poten- 
tials for A^A^ scattering, the most advanced potentials [H El [3l H] still use phenomeno logical 
form-factors to account for intermediate to short range interactions. Such form-factors shed 
little light on how nuclear forces at short distances follow from the basic concepts of QCD. 
Presently only the long range A^A^ interaction is understood on fundamental QCD grounds. 

The situation is not spectacular also from the experimental point of view. New experi- 
ments aimed at studies of A^A^ interaction at short distances practically stopped after the 
reassignment of AGS at Brookhaven National Laboratory [66j. 

In this respect an alternative way of studying nuclear forces at short distances is to 
probe A^A^ systems in nuclei at short space-time separations. Expectations that this can be 
achieved only at high- momentum transfer reactions (see e.g. [71 El El [101 [H]) was confirmed 
in a series of experiments with high energy electron[T2l [131 El lEl HE] and proton probes p!7l 
flSl flQl [20l [21] . Some of the unique results of these experiments were the observations 
of the scaling for the ratios of inclusive cross sections of nuclei and the deuteron [12j (or 
^iife[T3l [H]) at xsj > 1 at > 1.5 GeV^ as well as the observation of the strong (by 
factor of 20) dominance of pn relative to pp/nn short-range correlations in the ^^C nucleus 
for bound nucleon momenta 300-600 Me V / c [TBI . If the first result was an indication 
that two (or more) nucleons can be probed at small separations, the second one was an 
indication of the dominance of the tensor- forces [231 [21] in such correlations. 

The simplest reaction which could be used to investigate short-range A^A^ interactions 
using nuclear targets is the exclusive electrodisintegration of the deuteron in which large 
magnitudes of the relative momentum of the pn system in the ground state are probed. 
Three experiments [25] [261 [27] have already been performed using the relatively high (up to 
6 GeV) energy electron beam of the Jefferson Lab and more comprehensive experimental 
program will follow after the 11 GeV upgrade of the lab [28]. 

The prospect of having detailed experimental data on high energy deuteron electrodisin- 
tegration makes the development of theoretical approaches for description of these reactions 
a pressing issue. 

One of the first models for high energy break-up of the deuteron were developed in the 
mid 90's in which main emphasis was given to the studies of nucleon rescattering in the final 
state of the reactions[29, 301 [3ll [32]. These models were simple in a way that they assumed 
a factorization of the electromagnetic 7*A^ and final state A^A^ interactions and considered 
the rather small values of relative momenta of the bound pn system. 

The extension of these calculations to a larger internal momentum region required more 
elaborate approaches and several studies were done in this direction [TU] [55] [2S1 [2S1 [571 

In this work we develop a theoretical model for the description of high exclusive 
electrodisintegration of the deuteron in knock-out kinematics based on virtual nucleon ap- 
proximation. The main theoretical framework is based on the generalized eikonal approxima- 
tion (GEA)^ [311 [3511101 [m [121 [13] which allows us to represent the scattering amplitude 
in the covariant form using effective Feynman diagram rules. In this way all the virtualities 
involved in the scattering amplitudes are defined unambiguously. Reducing these amplitudes 
by choosing positive energy projections for the nucleon propagators allow us to represent 



2 



them through the convolutions of the deuteron wave function, on- and off- shell components 
of electromagnetic current and pn rescattering amplitude. In addition to accounting for 
the off-shell effects, nonfactorized approximation is applied to the electromagnetic and A^A^ 
rescattering parts in the calculation of the final state interaction (FSI) amplitude. As a 
result our calculation extends beyond the distorted wave impulse approximation limit. We 
also estimated the charge exchange contribution in the final state interaction in addition to 
the pn pn rescattering part of the FSI amplitude included in the eikonal approximation. 

The paper is organized as follows. In Sec. 2 we briefly discuss the kinematics of the 
disintegration reaction which we consider most efficient in probing the pn system at small 
separations. Then we discuss the basic assumptions of virtual nuclear approximation and 
proceed with the derivation of scattering amplitudes and the differential cross section of the 
reaction. 

In Sec. 3, after deriving the total scattering amphtude we performed a detailed theoretical 
analysis to identify the extent of uncertainties due to the off-shell part of the final state 
interaction as well as contribution due to charge-exchange rescattering. We also analyzed 
the role of the off-shell effects in the electromagnetic current of the bound nucleon. These 
analyses allowed us to conclude that at sufficiently large values of (~ 6 GeV^) the three 
most important contributions into the disintegration process are the off-shell electromagnetic 
current of the bound nucleon, the deuteron wave function and the on-shell part of the A^A^ 
scattering amplitude. 

Furthermore we compare our calculations with the first available high experimental 
data. These comparisons allow us to confirm the early prediction of GEA that the maximal 
strength of FSI corresponds to ~ 70° production of the recoil nucleon relative to the direction 
of the virtual photon. We also found that forward angles of recoil nucleon are best suited 
for studies of the off-shell electromagnetic current and the deuteron wave function. Another 
observation is that in backward direction there is a sizable contribution due to the A-Isobar 
production at the intermediate state of the reaction. In Sec. 4 we give conclusions and an 
outlook on further development of the model. 

II. CROSS SECTION OF THE REACTION 
A. Kinematics 

We discuss the process: 

e + d^e'+p + n, (1) 

in knock-out kinematics in which case one nucleon (for definiteness we consider it to be 
a proton) absorbs the virtual photon and carries almost all its momentum. The optimal 
kinematics for probing the initial pn system at close distances is defined as follows: 

(a) Q2 > 1 GeV^; (b) pf ^ q- (c) Pf :$> Pr > 300 MeV/c, (2) 

where we define q = {qo,q), Pf = {Ef,pf) and Pr = {Er,Pr) as four-momenta of virtual 
photon, knock-out proton and recoil neutron respectively. Also = — Q'o- Conditions 
(|2|(b) and (c) define the knock-out process, while condition ([2]) (a) is necessary to satisfy 
Eq.([2])(c). From the point of view of the dynamics of the reaction one also needs Eq.([2])(a) 
in order to provide a necessary resolution for probing the deuteron at small inter-nucleon 
distances. 
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FIG. 1: GEA diagrams 

In the most simple picture the kinematics of Eq.(|2]) represent a scenario in which the high 
energy virtual photon removes the proton from the pn system leaving the neutron with the 
pre-existing relative momentum pr (see Fig.([T])a). 

B. Main Assumptions of Virtual Nucleon Approximation 

The virtual nucleon approximation is based on the following main assumptions which also 
define the limits of its validity. First, one considers only the pn component of the deuteron, 
neglecting inelastic initial state transitions. Since the deuteron is in a isosinglet state this 
will correspond to restricting the kinetic energy of recoil nucleon to 

Tn < 2(mA - ttin) ~ {rriN* - m) ~ 500 MeV (3) 

where m, rriA and rriN* are masses of the nucleon and low-lying non-strange baryonic reso- 
nances, A(1232) and A^*(1525). We neglect also the pionic degrees of freedom in the deuteron 
wave function. However we expect that the overall error introduced by this approximation to 
be small since the probability of low momentum pionic degrees of freedom is suppressed due 
to pseudogoldstone nature of pions in QCD as well as observation that nNN form factors 
are soft (see e.g. discussion in Refs.[5| fTT]). 

The second assumption is that the negative energy projection of the virtual nucleon 
propagator gives negligible contribution to the scattering amplitude. Such an assumption 
can be justified if, 

Md - ^jw? + p2 > 0, (4) 

where is the mass of the deuteron and p is the relative momentum of the bound pn 
system. 

The above two conditions can be satisfied if we restrict the momentum of the recoil 
neutron, p^. < 700 MeV/c. However due to the fact that we explicitly left out the non- 
nucleonic components of deuteron wave function, the momentum sum rule is not satisfied 
in virtual nucleon approximation (see discussions in Refs.([8l 1591 160] ). 

The third assumption which is made in the calculations is that at large (> 1 GeV^) 
the interaction of virtual photon with exchanged mesons are a small correction and can be 
neglected (see e.g. discussions in Ref . P [TO] ) . 
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C. Generalized Eikonal Approximation 



The assumptions discussed above allow us to restrict the consideration by the set of 
Feynman diagrams presented in Fig.Q. One can calculate these diagrams based on the 
effective Feynman diagram rules discussed in Ref.[Tn]. These rules allow us to formulate 
scattering amplitudes in the covariant form which unambiguously accounts for all the off-shell 
effects. Then we reduce the covariant amplitudes into a non-covariant form by choosing the 
positive energy projection of nucleon (or baryonic resonance) propagators at the intermediate 
state of the reaction. 

Figjlj^a) diagram corresponds to the plane wave impulse approximation (PWIA) in which 
the virtual photon knocks out one of the nucleons from the deuteron leaving the second nu- 
cleon in the on-shell positive energy state. Two nucleons do not interact in the final state 
of the reaction representing two outgoing plane waves. The diagram of Fig|l](b) represents 
a situation in which the elastic electroproduction is followed by the elastic pn — > pn rescat- 
tering. In this case the rescattering is forward in the sense that, for example, the proton 
struck by the virtual photon will attain its large momentum after the pn — »• pn rescattering. 
The amplitude of this scattering will be referred to as a forward FSI amplitude. 

The diagram of Fi^TJ^c) corresponds to the scenario in which final state interaction pro- 
ceeds through the charge-exchange pn — > np rescattering. In this case the final fast proton 
emerges from the process in which the photon strikes the neutron which then undergoes 
a. np ^ pn charge-exchange rescattering. The amplitude of this scattering process will be 
referred to as a charge-exchange FSI amplitude. 

The fourth diagram (FigjlJ^d)) corresponds to the electroproduction of an excited state 
R with a subsequent RN NN final state rescattering. The most important contribution 
to the fourth diagram is due to the A-isobar (IC), whose production threshold is closest to 
the quasielastic scattering kinematics. Several factors make IC contribution small in high 

limit at X > l[9l [10]. One factor is the large longitudinal momenta of the initial nucleon 
involved in the A-electroproduction process: 



2 2 

= (1 - x)m - , (5) 

which indicates that choosing x > 1 one can suppress the electroproduction of A-resonance 
in the intermediate states due to the large values of initial momenta entering in the deuteron 
wave function. 

An additional suppression of IC follows from the smallness of the 7*iV — > A transition 
form-factors as compared to the elastic form factors at > few GeV^[l9t 150] . Finally, 
due to the fact that the AA^ NN amphtude is dominated by pion or p-type reggeon 
exchanges, it will be additionally suppressed by at least the factor of — p=|- In any case 

this contribution can be calculated in a selfconsistent way within the generalized eikonal 
approximation. The calculation of these types of diagrams within GEA will be presented 
elsewhere [HJ. 

Below we will discuss the calculations of only the first three diagrams of Figjlj 
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1. Plane Wave Impulse Approximation Amplitude 

The amplitude of the PWIA diagram in the covariant form can be written as follows: 

{Sf, Sr\A^\ Sd) = -u{Pr, Sr)T'^,p ^^^^^ ■ u{pf, Sf)rDNN " (6) 

Pi 

where r^*p is the electromagnetic vertex of the 'y*N ^ N scattering and the vertex function 
^DNN describes the transition of the deuteron into the pn system. The notations Sd,Sf,Sr 
describe the spin projections of the deuteron, knock-out proton and recoil neutron respec- 
tively. The spin function of the deuteron is represented by x^'*- The four-momentum of the 
struck nucleon in the initial state within PWIA is defined as: 

Pi = {Ed - Er,Pd - Pr) = {Md - Er, -pr) \LaB ■ (7) 

The above relation clearly shows the off-shell character of the struck nucleon in the initial 
state, since pj 7^ m^. Therefore expressing the initial nucleon's propagator through the 
on-mass shell nucleonic spinors is not valid. 

However, using an approximation in which only positive energy projections are taken 
into account, one can isolate the on-shell part of the propagator by adding and subtracting 

J^on^O |;gj~j-Q f^Q r^. g^g follows: 

^ + m = + m + {Ef^ - E°")7°, (8) 



where E°^-^ = Md — \J + p^ and E""^ = ^vriy + pi where m„ and are the masses of the 
proton and neutron respectively [67]. Now we can separate the PWIA amplitude into on- 
and off-shell parts in the following way: 



(S/, Sr\Al\ Sd) = {Sf, Sr \ A^,^ \ Sd) + {Sf, Sr \ A^ \ Sd), (9) 



where 



and 



{Sf, Sr I v4^,o„ I Sd) = -U{pr, Sr)T^*p ^l J^" ^ " u{pf, Sf)TnNNXd, (10) 

Pi ^ 

{Sf, Sr I I Sd) = -U{pr, Sr)T^,p ' 2 _ \ ^^^Z' ^f)^DNNXd- (H) 

Pi ^ 

For the on-shell part of the amplitude, using 

+ m = J2u{Pi,Si)u{pi,Si) (12) 

Si 

and the definition |^ H6] . 

nis.P.S2,p2) = - ^ (13) 

{pj - m2)v^^(27r)3(p2 + ^2)1 

one obtains 



Sf,Sr I ^0,o„ I Sd) = V2j{2Txf2ErY,JN,on{Sf,Pf]SuPi)%\Si,Pi,Sr,Pr), (14) 
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where 

JN,on{Sf,Pf;Si,Pi) =u{pf,Sf)T^*^u{pi,Si). (15) 

For the off-shell part of the scattering amplitude one observes that the following relation, 
u{p2,S2)TdnnX = -2^ 7. {Pi -"^ )v2V(27r)^2(pf + m^)2 (16) 



satisfies Eq.(13). Inserting it into Eq.(ll) one obtains 



{Sf,Sr I A^^^fj I Sd) = V2^ {2Tlf2Er ■JN,offiSf,Pf; Si,Pi)'^'/{Si,Pi, Sr,Pr), (17) 

Si 

where 

2m 



JN,off{Sf^Pr^Si^Pi) =u{Pf,Sf)T^*Nl°u{pi,Si)^— (18) 



and Ef^ = Md - E^ and = ^m^ + pi 

One can combine on-and off-shell parts of the PWIA scattering amplitudes in the following 
form: 



{Sf,Sr I I Sd) = V2^{27ry2ErY,J^iSf,Pf;Si,Pi)'^'/{Si,Pi,Sr,Pr), (19) 

Si 

where 

JNisf,Pf,Si,Pi) = J^^^^{sf,pf;Si,Pi) + J^^^ff{sf,pf;Si,Pi). (20) 



The above form of the electromagnetic current together with Eqs.(15) and ([18j) represents 
our off-shell approximation. It is worth noting that the first, "on-shell" part of this current 
corresponds to the widely used "De Forest" approximation [17] . In the "DeForest" approxi- 
mation because of the absence of the second term the gauge invariance is violated and the 
current conservation is restored by expressing Jq or components through each other with 
different assumptions for the nucleon spinors and electromagnetic vertices. The latter intro- 
duces more uncertainty since imposed relations are not unique. As a result one generates 
several forms of the off-shell electromagnetic currents |4 7] . 



The additional "off-shell" part of the electromagnetic current in Eq.(20) obtained in 



our approximation reduces the uncertainty of choosing on-shell nucleon spinors which is 



inherent to the "De Forest" approximation. The total current in Eq.(20) is conserved since 
it is derived from the gauge invariant amplitude of Eq.([6]). Therefore our approximation 
does not violate gauge invariance and no additional conditions are needed to restore the 
current conservation. 

Note that our approximation is analogous to the one used in hadronic physics within 
light-cone approximation (see e.g. |18]) in which case an off-shell "7+" component of the 
fermion propagator is isolated in the similar manner as it is done for the 7° component in 
our case. 



2. Forward Elastic Final State Interaction Amplitude 

We start by applying effective Feynman diagram rules to the diagram of Fig.([T|)b, which 
yields: 



{Sf,Sr \ A'^ \ Sd) = - J 



d P'r U{pf, Sf)u{pr, S,)FNNWr + m][pd - Pr + + 

i{2T{Y [pd - p'r + qY -rn? + ie 



m\ 
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{{Pd — p'r^ — m? + 'i^){p? — rn? + ie) ' 

were Fnn represents the invariant pn — > pn scattering amplitude that can be expressed as 
follows: 

Fnn{s, t) = ^s{s-Am^)fNN{s, t) (22) 

where s is the total invariant energy of the scattering pn system and the /at at scattering 
amplitude is defined in such a way that Imfj^jsi = atot- Furthermore we will use the following 
four-vectors defined as 

P'i =Pd- p'r and p'f =Pd- p'r + q- (23) 

We first integrate by d^p'^ through the positive energy pole of the spectator nucleon propa- 
gator at the intermediate state; 

This integration allows us to use "/ir + m = J2 u{p'^, s'^)u{p'^, s',). For — f we use a relation 



similar to Eq.(|8]). The same could be done for "pa — p[. + ^. However for large values of 
q the off-shell part in Eq.(8) is suppressed by ^^^j^ and in large limit it's contribution 
is negligible. Thus we can use the on-shell relation, i)d — i>' ^ k — J2u{p'f, s'j-)u{p'j, s'jr) for 

the knock-out nucleon spinors in the intermediate state. Using the above representations 
of the spinors, the definitions of the deuteron wave function (Eq.(13)) and electromagnetic 
current (Eq.(20)) for the scattering amplitude of Eq.(21) we obtain: 



3 r dV j2E'.Js(s -Am'^) 

J i{27ry 2E'^{{pd -p'r + qr -m^ + le) 



{pf,Sf]Pr,Sr I / {t,s) I p^,s'^]p'f,s'j) ■ J'^{s'f, p'f] Si, Pi) ■ %\Si, p[, s[., p[) . (25) 

Next, we consider the propagator of the knock-out proton in the intermediate state, using 
the condition of quasielastic scattering 



one obtains 
where 



{q+Pd-Pr) =Pf = 'm, (26) 



{Pd-p'r + q) + = 2|q|(p^^2 -pr,2 + A + ze), (27) 

A = ^(E,.-E-M!^(E,.-E;.) + '!t^. (28) 



Furthermore using the relation 
1 
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and performing integration over ^ we split the scattering amplitude into two terms; one 
containing on-shell and the other off-shell pn pn scattering amplitudes as follows: 



' ' ' ' ' ' ^ " ^^4^./ ^^^^^ 

{pf,Sf;Pr,Sr I I ■ ■ ^?(^»'^i'4>Pr) 

X 



2 7 27r J (27r)2 2E;|q 



(p/, 3/; p„ Sr I f^^'"ff{t, s) I p;, g;; p}, 4) , , , , 
, _ —JN\Sf,Pf]SuPi) ■ "^/{Si^p^^s^^p^}, (30) 

Pr,z Pr,z 



where p'^ = {pr,z - A,p;_l), E'^ = ^m? + p'^, p'- =p^-p'^ and p'^ = p'^ + q. 

For numerical estimates of the above amplitudes one needs on- and off-shell pn pn 
amplitudes as an input. In high energy limit in which the helicity conservation of small angle 
NN scattering is rather well established the on-shell amplitude is predominantly imaginary 
and can be parameterized in the form 

(p/, Sf\pr, Sr I /^^'°"(t, s) I p^, s[.]Pf, s'f) = + a)e^^6sf,s'/sr,s'^, (31) 

where crf^(s), B{s) and a{s) are found from fitting of experimental data on elastic 
pn — > pn scattering. For the effective lab momentum range of up to 1.3 GeV/c the SAID 
parameterization [5 Ij of pn amplitudes can be used. The situation is more uncertain for the 
half-off-shell part of the f^^'°fi amplitude. In present calculations we use the following 
parameterization: 

fNN,off ^ jNN,on^B{nilj^-m-^) ^ ^gg) 

where '^o// = — E'^ + qoY — (p^ -|- g)^. Overall we expect that our calculation will not 
be reliable in situations in which the contribution from the off-shell part of the rescattering 
is dominant. However in high limit this contribution is only a small correction. 

Completing this section it is worth to notice that in addition to the appearance of the 



A factor (Eq.(28)) in GEA which does not enter in conventional Glauber approximation 



(see detailed discussion in Ref. [10]), the new factor, '^^2E'\q\' " entering the elastic FSI am- 



plitude (Eq.(p30|) is also unique to GEA. Within conventional Glauber approximation, in 
which Fermi motion of the scatterers is neglected this factor is equal to one. However within 
GEA it appears as a consequence of the covariant form of the initial scattering amplitude. 
Calculation of this factors for our kinematics yields: 



sis-Am^) ^C-^Q^-mDC-^ ^^^^ 



2E'^\q\ 2E'^\q\ 

which decreases with a; — >■ 2. Thus for the a; > 1 and large kinematics GEA predicts an 
additional suppression of FSI as compared to the conventional Glauber approximation. 
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3. Charge- Exchange Final State Interaction Amplitude 



To complete the calculation of the total amplitude we need to include the contribution 
from charge-exchange rescattering, which can be obtained from Eq.(30) after the substitu- 
tions corresponding to Fig.lc. Namely, one needs to switch the proton and neutron lines in 
the initial and intermediate states of the scattering, replace proton electromagnetic current 
by the neutron and fNN by the charge-exchange scattering amplitude fj^^ ■ One obtains: 



j\/2(27r) 



E 



pE'Js{s-Am^) 



(27r)2 2E',\q\ 



27T 



(2vr) 



2E'M\ 



{pf,Sf,Pr,Sr I \p'^,s'/,p'f,s'f) 
Pr.z Pr.z 



Ji:{s'j,p'f;s„p[)-^l^{s„p[,s',,p',). {34) 



Here the charge-exchange rescattering amplitude similar to the elastic FSI case is taken 
from the experimental measurements. The off-shell extrapolation of the rescattering ampli- 
tude is also done similar to Eq.(32). For numerical estimates we use the parameterization 
of Ref.l 



4- Deuteron Wave Function 



The deuteron wave function in Eq.(13) in general represents a solution of the Bethe- 



Salpeter type equation. To fix the normalization of the wave function we need to relate 



an expression that contains the deuteron wave function (as it is defined in Eq.(13)) to a 
well defined quantity characterizing the deuteron. One such quantity is the deuteron elastic 
charge form-factor Gc, which at — >• limit approaches to one, i.e. Gc{Q^ = 0) = 1 (see 
e.g. Ref. [53]). The latter could be related to the deutron elastic scattering amplitude as 
follows: 

^ jZ {p'„ s', I A^=\Q') I p,, |q2_o= Gc(0) = 1, (35) 

where {p'd,s'^ \ \ pa, Sd) is the elastic 7*c? d' scattering amplitude corresponding to the 
diagram of Figj2] 

Applying the same effective Feynman diagram rules used above for (p'^, s'^ \ \ pd, Sd) 
one obtains: 

J4 



/ / / I I \ ^ f d'^Pr s' trt h + ^ -rf. i>i + m 



X . . . 36 

pf. — + le 

Further derivations within the virtual nucleon approximation follow the similar to 
Secs.II.C.l and II. C. 2 steps. We first evaluate integral by the pole value of the spectator 
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FIG. 2: Elastic 'y*d —>■ d! diagram. 



nucleoli propagator, then separate the on- and off-shell parts in the numerator of interacting 



nucleon propagators and introduce deuteron wave function according to Eq.(13). In this 



case the electromagnetic current of the nucleon is fully off-shell since the struck nucleon is 
bound in both initial and final states of the reaction. This results in the following expression 
for the elastic scattering amplitude: 



El^f - El 



1 + 



- ^2°" 

2m 



-70 



J- ^*Ar 



/ + 



2m 



-70 



u{pi,Si)'$'/{Si,Pi,Sr,Pr). (37) 



Neglecting the second order off-shell terms in the above equation (i.e. 
obtains: 



' E°ff-E°" \2^ 
2m > > 



l,n S2,Sl,Sr 



{P'd^s'dl A^" \ Pd,Sd) = 4:J2 E I d^Pr'^'f\s2,P2,Sr,Pr)JN'^'d{Si,Pi,Sr,P 



one 



(38) 



where 



Jn{s2,P2] Si, Pi) = {S2,P2; SuPi) + JN,off{s2,P°2^] Si,pi) + J^„ff{s2,P2; Supf^). (39) 



Here the on- and off- shell parts of electromagnetic current are defined in Eq.(15) and (18). 
In the above equation the notation p"'^-!' in the argument of JN,off indicates which nucleon 
is considered as off-shell. 

Using now the fact that for the proton and neutron -Fip(„)(Q^ = 0) = 1(0) and using 
Eqs.( p|l8| one obtains: 



=0 



2El^^ 
0. 



Using these relations and inserting Eq.(38) into Eq.(35) one obtains 

1 



(40) 



(41) 



where E°f^ = Md — \Jm^ + p^. It is worth mentioning that above normalization coincides 
with the normalization obtained from the baryon number conservation sum rule[5ll |55l [56l 
[571 158| [59| 160] for deep inelastic scattering off the deuteron: 



(42) 
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where a = is the hght-cone momentum fraction of the deuteron carried by 

the struck nucleon. As it was mentioned in Sec. II. B in virtual nucleon approximation due 
to unaccounted non-nucleonic degrees of freedom the wave function defined according to 
Eq.(13) will not satisfy the energy-momentum sum rule which expresses the requirement 
that the sum of the light-cone momentum fractions of all the constituents of the nucleus 
equals to one. 

For numerical estimates we model the deuteron wave function to satisfy Eq.(41) in the 
following form [561 160] : 

^dip) = ^r(p) ,,,. % Tr^ (43) 

2{Md - vm^ +p^) 

which provides a smooth transition to the nonrelativistic wave function "^^^ in the small 
momentum limit. 



5. Total amplitude and the differential cross section 

The total scattering amplitude consists of the sum of PWIA, forward and charge-exchange 
FSI amplitudes: 

(S/,Sr I A^" I Sd) = {Sf,Sr I A'^ I Sd) + {Sf,Sr \ A^ \ Sd) + {Sf,Sr \ I ^d)- (44) 

Using this amplitude the differential cross section is calculated as follows: 

_a^K I I J^Jd,, ? P) , . 



where 

= u{k2,S2)^^u{h,Si) (46) 

and 



/2(27r)32E, 

For numerical estimates we use the electromagnetic current of the nucleon in the form 

F'^ = F,{Q')^^ + ^^^ra^^^q,, (48) 

where Fi and F2 are Dirac form-factors and for their evaluation the available phenomenologi- 
cal parameterizations are used[52]. For the deuteron wave function we use the approximation 



of Eq.(43) with the non-relativistic wave function calculated based on the Paris potential[l]. 
The pn scattering amplitude is parameterized in the form of Eq.(31 ) and its off-shell extrapo- 
lation in the form of Eq.(32). Also, for /p„ in the lower momentum range {piab ^ 1-3 GeV/c) 
we use the SAID parameterization [51] based on the pn scattering phase-shifts. 
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III. OBSERVABLES 



The main quantity which we will consider for numerical estimates is the ratio of the calcu- 
lated cross section which includes total amplitude of Eq.(44) to the cross section calculated 
within PWIA: 

R = -Z^A (49) 



where a == dn^idpidiif ■ This ratio allows us to clearly distinguish between kinematics 
in which PWIA dominates R 1 from kinematics in which FSI is dominated by screening 
i? < 1 or single rescattering R> 1 effects (see e.g. SHI E2])- 

Considering the numerical estimates of the ratio R we will discuss four main effects which 
characterize our present theoretical approach. These are the unfactorization of the electro- 
magnetic interaction in the FSI amplitude, the off-shell effects in the final state interaction, 
the effects of charge-exchange rescatterings and the off-shell effects in the electromagnetic 
interaction of the bound nucleon. 

In our estimates we will study the dependence of R on the angle of the recoil neutron 
relative to q for different values of neutron momenta. We will perform our calculations for 
two values of (2 GeV^ and 6 GeV^) which will allow us to also assess the dependence 
of the considered effects. 

Finally, we will present the comparisons with the first experimental data on the deuteron 
electrodisintegration at large Q^. 




FIG. 3: Dependence of ratio R on the recoil angle of the neutron for different values of pr = 
100,200,300,400,500 MeV/c and = 2,6 GeV^. Solid line - unfactorized and dashed line - 
factorized approximations. 
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A. Nonfactorization effects 



In Fig j3] we compare the calculations of R with and without factorization approximation 
for the electromagnetic current in the FSI amplitude. The factorization approximation will 
result in the so called distorted wave impulse approximation (DWIA) widely used in the 
literature. 

As the figure shows the factorization (or DWIA) approximation is applicable for up to 
Pr < 300 MeV/c or for the kinematics in which the FSI amplitude is smaller than the 
PWIA term. The factorization approximation breaks down in kinematics dominated by 
the rescattering process at pr > 400 MeV/c. As the comparisons show, the unfactorization 
predicts a larger FSI amplitude which can be understood based on the fact that in this case 
the electromagnetic current which enters in the rescattering amplitude of Eq.(30) is defined 
at smaller values of bound nucleon momenta than the electromagnetic current in the PWIA 
term (Eq.(|l9[). 

Figjs] shows also that the factorization approximation is dependent and somewhat 
improves with an increase of Q^. This is a rather important feature which should be taken 
into account in color transparency studies (CT) for double scattering kinematics, in which 
case the dependence of the peak of the ratio R is studied in order to extract the CT signal 
(see e.g [311 162|, [63]). Our comparisons in Figj3] shows that the unfactorized approximation 
should be used as a baseline for identification of the CT signature in the dependence of 
the FSI contribution of the deuteron break-up cross section. 




FIG. 4: Dependence of ratio R on the recoil angle of the neutron for different values of pr and 
Q2 = 2,6 GeV2. The recoil neutron momenta are the same as in Fig. 3. Solid line - unfactorized 
calculation with the pole term only in the FSI amplitude, dashed line - unfactorized approximations 
with both pole and principal value terms in the FSI amplitude. 
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FIG. 5: Dependence of ratio R on the recoil angle of the neutron for different values of pr and 
= 2,6 GeV^. The recoil neutron momenta are the same as in Fig. 3. Solid line - unfactorized 
approximation with the pole term only forward pn pn rescattering, dashed line - unfactorized 
approximation including the pole terms for both forward pn pn and charge-exchange pn — > np 
rescattering amplitudes. 



B. Off-Shell FSI Effects 

Next we consider the contribution to the FSI amplitude due to the principal value part of 



Eq.(30). This term depends on the half-off shell A^A^ scattering amplitude which is a largely 
unknown quantity. Therefore the reliability of our calculations depends on the magnitude 
of the principal value term. In Fig|4] we compare the calculations in which only the pole 
(on-shell) term of the FSI amplitude is included with the calculations which contain both 
pole (on-shell) and principal value (off-shell) terms of the FSI amplitude. For the half-off 



shell fpn amplitude we use the approximation of Eq.(32). An important observation can be 
made from Fig|4] is that the off-shell rescattering effects diminish with an increase of Q^. 
This is consistent with our earlier observation |1 11 [12] that the distances that the virtual 
nucleon propagates before the rescattering decreases significantly with an increase of at 
fixed values of x. 



C. Charge Exchange Rescattering Effects 



Due to the fact that the charge-exchange rescattering amplitude is predominantly real it 
will interfere mainly with the real part of the forward elastic FSI amplitude which is a small 
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FIG. 6: Ratio of PWIA cross sections calculated using on-shell and off-shell approximations for 
the electromagnetic current of the proton. The curves from top to bottom correspond to recoil 
neutron momenta 100, 200, 300, 400, 500, 600 and 700 MeV/c respectively. 

parameter at large energies. One can see from Fig|5]that the charge-exchange rescattering 
dominates at kinematics in which the rescattering defines the overall magnitude of the cross 
section. 

However it is a rather well known fact that, due to the dominant pion-exchange nature 
of charge-exchange pn scattering, its cross section decreases linearly with an increase of the 
invariant energy s as compared to the forward pn elastic scattering cross section. As a result 
one expects that with an increase of the charge exchange rescattering term will become 
a small correction. This can be seen from the calculation for = Q GeV^ kinematics in 
FigjSj 

D. OfF-Shell Electromagnetic Interaction Effects 

The above evaluations of the d{e, e'p)n cross sections demonstrate that in the large 
limit the property of the scattering is defined mainly by the PWIA and forward angle on- 
shell FSI terms. The angular distribution has a very characteristic shape in which one can 
identify kinematics dominated by PWIA or FSI. Calculations also show that the forward or 
backward kinematics of the recoil nucleon are best suited for isolating PWIA term which 
subsequently allows us to study the structure of the electromagnetic current of the bound 
nucleon and the deuteron wave function at large values of internal momenta. 

We now concentrate on the effects related to the fact that the proton in the deuteron 
becomes increasingly off-shell at large values of recoil neutron momenta in forward/backward 



= 2 GeV^ 
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directions. 

As it follows from Eqs.( l5p!8 ) and (20) the off-shell part of the electromagnetic current 
will diminish the overall magnitude of the electromagnetic interaction. Since the off-shellness 
grows with an increase of the initial momentum of the struck nucleon it will result in the 
suppression of the electromagnetic interaction strength of deeply bound nucleons. As it 
follows from Fig|6]the off-shell effects have a weak dependences and to disentangle them 
from the effects related to the high momentum component of the deuteron wave function 
would require measurements covering an extended angular and recoil momentum range. 
Figj6] also shows that the forward direction of recoil nucleon momenta represents the most 
optimal kinematic region for minimizing the off-shell effects for electromagnetic interaction. 

Note that polarization measurements will provide additional observables for separating 
current and wave function effects. For example, measurements of the cross section asymme- 
tries similar to Ref. at large recoil momenta will be more sensitive to the structure of the 
electromagnetic current since wave function effects largely cancel out in the ratios defining 
the asymmetry. 



E. Comparison with Experimental Data 

In the last few years three experiments P5| f2E\ [27] produced first data at relatively large 
kinematics. 

The first experiment [25] probed the = 0.665 GeV^ and x ^ 1 kinematics and provided 
very accurate data. The measured value of is marginal for the application of GEA. 
However as Fig j7] shows the comparison yields a surprisingly good agreement with the data. 
FigjT] compares the reduced cross section defined as follows: 



where the differential cross section is defined according to Eq.(45) and acci is the off- 
shell electron-proton cross section defined in Ref. |47j . Note that in this calculation for 
fpn amplitude we use SAID[51j parameterization for both elastic and charge-exchange pn 
scatterings which fit the elastic pn scattering data for lab momenta up to 1.3 GeV/ c. Because 
of the relatively low energy and momentum transfers involved in the reaction the calculation 
shows a substantial contribution from the off-shell as well as charge exchange parts of the 
FSI amplitude (difference between dotted, dash-dotted and solid lines in Figj?]). 

It is worth noting that in agreement with the observations of Refs.[27| 138] the calculation 
overestimates the low recoil momentum part of the cross section (FigjTjD), where we expect 
that theoretical uncertainties are well under control. Our preliminary estimates demonstrate 
that this discrepancy could be accounted for by inclusion of the contribution of two-photon 
exchange effects in the overall amplitude of the scattering |65j. 

The second experiment |26]) covered the range from 2 — 5 GeV^. However, due to the 
low statistics the data are integrated over the ranges of the recoil nucleon momenta. 

Comparing with these data we performed the same kinematic integrations as the experi- 
ment did. The results are shown in Fig|8]and[9j Despite these integrations we still can make 
several important observations from these comparisons. 

• First, the angular distribution clearly exhibits an eikonal feature, with the minimum 
(Figjsl) or maximum (Figjo]) at transverse kinematics due to the final state interaction. 
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FIG. 7: Missing momentum dependence of the reduced cross section. The data are from Ref . |25] . 
Dashed Hne - PWIA calculation, dotted line - PWIA+ only pole term of forward FSI, dash-dotted 
line - PWIA-I- forward FSI, solid line - PWIA + forward and charge exchange FSI, and solid line 
with squares - same as the previous solid line, added the contribution from the mechanism in which 
the proton is a spectator and the neutron was struck by the virtual photon. 

The most important result is that the maximum of FSI is at recoil angles of 70'' in 
agreement with the GEA prediction of Ref. |10] . Note that the conventional Glauber 
theory predicted 90° for the FSI maximum. 

• The disagreement of the calculation with the data at Or > 70° appears to be due to the 
isobar contribution at the intermediate state of the reaction. This region corresponds 
to X < 1 and it is kinematically closer to the threshold of A-isobar electroproduction. 
The comparisons also indicate that the relative strength of the A-isobar contribution 
diminishes with an increase of and at neutron production angles 6r — > 180°. 

• The forward direction of the recoil nucleon momentum, being far from the A-isobar 
threshold, exhibits a relatively small contribution due to FSI. This indicates that the 
forward recoil angle region is best suited for studies of PWIA properties of the reaction 
such as the off-shell electromagnetic current and deuteron wave function. 

Finally, the results of the experiment of Ref. [27] are currently in the final stages of analy- 
sis. Since in this experiment the disintegration reaction is measured at forward recoil angles 
and at up to 3.5 GeV^, we expect that it will allow us to further check the validity of 
our claim that the forward angular region is best suited for studies of the PWIA properties 
of the reaction. 
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FIG. 8: Dependence of the differential cross section on the direction of the recoil neutron momen- 
tum. The data are from Ref . |26j . Dashed line - PWIA calculation, dotted line - PWIA-I- pole term 
of forward FSI, dash-dotted line - PWIA-l-forward FSI, solid line - PWIA -|- forward and charge 
exchange FSI. The momentum of the recoil neutron is restricted to 200 < pr < 300 MeV/c. The 
labels 2, 3, 4 and 5 correspond to the following values of = 2±0.25; 3±0.5; 4±0.5; 5±0.5 GeV^. 

IV. CONCLUSIONS 

Within the virtual nucleoli approximation we developed a theoretical framework for cal- 
culation of high Q"^ exclusive electrodisintegration of the deuteron at large values of recoil 
nucleon momenta. The scattering amplitude is derived based on generalized eikonal approx- 
imation, in which case each amplitude is calculated based on effective Feynman diagram 
rules. Because of the covariant formulation of the problem the electromagnetic current is 
gauge invariant from the beginning. By isolating the off-shell part in the electromagnetic 
current we introduced an approach which allows us to express the bound nucleon current 
separately through the on-shell and off shell parts. 

Next, we derived the final state interaction amplitude based on GEA and expressed 
it through the on-shell and off-shell rescattering parts. No factorization approximation is 
assumed in the calculation of the FSI amplitude. The calculation of FSI also includes the 
amplitude due to the charge-exchange final state interaction. 

We performed numerical analysis of the obtained formulae to identify the level of un- 
certainties due to the factors included in the calculations. Overall, our conclusion is that 
with an increase of all the uncertainties related to the off-shell FSI and charge exchange 
rescattering became small and the total scattering amplitude is determined by the PWIA 
and on-shell elastic A^A^ rescattering. 
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FIG. 9: Dependence of the differential cross section on the direction of the recoil neutron momen- 
tum. The data are from Ref . |26j . Dashed line - PWIA calculation, dotted line - PWIA-I- pole term 
of forward FSI, dash-dotted line - PWIA-l-forward FSI, solid line - PWIA -|- forward and charge 
exchange FSI. The momentum of the recoil neutron is restricted to 400 < pr < 600 MeV/c. The 
labels 2, 3, 4 and 5 correspond to the following values of = 2±0.25; 3±0.5; 4±0.5; 5±0.5 GeV^. 
The data sets and calculations for "4" and "5" are multiplied by 0.5 and 0.25 respectively. 

We compared our calculations with the first experimental data at large deuteron 
electrodisintegration. These comparisons revealed a rather surprising agreement with low 

= 0.665 GeV^ data which indicates the wider range of applicability of the present ap- 
proximation. 

Comparisons with higher data (> 2 GeV^) at the wider range of recoil nucleon mo- 
menta and angles demonstrate the important role that the intermediate A- Isobar production 
plays in electrodisintegration reaction at backward angles close to the A production thresh- 
old. 

However, at forward recoil angles where FSI effects are restricted, the calculations show 
greater sensitivity to the PWIA structure of the electrodisintegration reaction. Further 
experiments will allow us to confirm that this region is most suitable for probing the bound 
nucleon electromagnetic current and the deuteron wave function at small distances. 
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